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functions involving the stress-energy tensor and other conserved currents. We then apply the 
obtained results on the correlation functions to computing the energy-energy correlations, which 
measure the flow of energy in the final states created from the vacuum by a source. We demon¬ 
strate that they are given by a universal function independent of the choice of the source. Our 
analysis relies only on A/" = 4 superconformal symmetry and does not use the dynamics of the 
theory. 
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1 Introduction 

The correlation functions of stress-energy tensors are very natural quantities to study in any 
conformal field theory. It is well known that the two-point function is fixed by conformal symmetry 
whereas the higher point correlation functions can have a very complicated form. Namely, they 
are given by a linear combination of numerous Lorentz structures arising as solutions to the 
conformal Ward identities consistent with the conservation of the stress-energy tensor. Their total 
number depends both on the number of points and the dimension of the space-time. In particular, 
in a four-dimensional conformal field theory, the three-point correlation function involves three 
different Lorentz structures [1-3]. They coincide with the three-point correlation functions of 
stress-energy tensors in a free theory of scalars, fermions and gauge fields, respectively. Eor 
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four-point functions, the situation is much more complicated - not only the number of Lorentz 
structures increases to 22, but each of them involves some function of the conformal cross-ratios 
[4]. This makes the calculation of the four-point correlation function of stress-energy tensors an 
extremely difficult task. To the best of our knowledge, there exists no closed expression for such 
correlation function in the literature. 

The problem becomes significantly simpler if the underlying conformal theory has super- 
symmetry. In general, it leads to additional constraints on the correlation functions of the 
stress-energy tensors and, therefore, greatly reduces the number of independent functions of cross¬ 
ratios. ^ In the four-dimensional maximally supersymmetric AA = 4 theory this number shrinks to 
one, so that the four-point correlation of stress-energy tensors depends on a single scalar function. 
The same function determines the four-point correlation function of the scalar 1/2 BPS operators 
020'- The reason for this is that the two operators, O 20 ' and belong to the same Af = A 
supermultiplet, the so-called stress-energy supermultiplet T, and their correlation functions are 
related to each other by = 4 superconformal Ward identities. The general solution to these 
identities, defining the four-point correlation function of the supercurrents (T(l).. .T(4)), was 
derived in [11]. Different correlation functions involving O 20 ' and appear as its components. 

In this paper, we derive the explicit expression for the four-point correlation function of 
stress-energy tensors in an AA = 4 superconformal theory. We show that it has a remarkably 
simple and suggestive form allowing us to extend the obtained results to a larger class of four- 
point correlation functions involving the stress-energy tensor and other conserved currents. Our 
analysis relies only on Af = A superconformal symmetry and does not use the dynamics of the 
theory. In the special case of A/" = 4 supersymmetric Yang-Mills theory (SYM), the four-point 
correlation functions are determined by the scalar function <h(tt,u), which has been extensively 
studied at weak and at strong coupling. 

Among all four-point correlation functions described by (T(1)T(2)T(3)T(4)), those involving 
two stress-energy tensors play a special role. They can be used to compute interesting (infrared 
safe) observables, the so-called energy-energy correlations [12], measuring the flow of energy in 
the final states created from the vacuum by some source J{x) [13-17]. The choice of the source 
is arbitrary and physically interesting cases involve the 1/2 BPS operator O 20 ', the stress-energy 
tensor, the ii—symmetry current and the Lagrangian of the theory, all belonging to the stress- 
energy supermultiplet T. The last two cases are of particular interest since they can be thought 
of as prototypes of the electromagnetic current in QCD and of the effective coupling of Higgs 
boson to gluons in the Standard Model, respectively. 

The energy-energy correlations can be obtained from the four-point correlation function 
(0|T^^i,^(l)r^2i"2(2)T(3)J(4)|0) through a limiting procedure described in great detail in [18, 19]. 
Using this approach, the energy-energy correlations have been computed \n Af = A SYM theory 
for a scalar source J = O 20 ' at next-to-leading order both at weak and at strong coupling [20-22] . 
Based on general physical considerations, we would expect that the flow of energy in the final state 
should depend on the choice of the source. Using the obtained results on the correlation functions, 
we find that, quite surprisingly, the energy-energy correlations \n Af = A superconformal theory 
are the same for the different source operators mentioned above and, therefore, are universal. 

The paper is organized as follows. In section 2 we describe the properties of the stress-energy 
supermultiplet and present the general expression for the four-point correlation function satisfying 

^For some general results on superconformal correlation functions see, e.g., [5-10]. 
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the AA = 4 super conformal Ward identities. In section 3 we explain how to extract the four-point 
correlation function of stress-energy tensors from the supercorrelator (T(1)T(2)T(3)T(4)). Then 
we use the special properties of the A/” = 4 superconformal generators to derive its explicit 
expression. In section 4 we extend the analysis to correlation functions involving conserved 
currents. We argue that they have a remarkably simple universal form and illustrate this by a 
few examples. In section 5 we apply the obtained expressions for the correlation functions to 
the evaluation of the energy-energy correlations for different source operators. Section 6 contains 
concluding remarks. Two appendices contain technical details. 


2 Superconformal Ward identities 

In this section, we summarize the properties of the correlation functions of the AA = 4 stress- 
energy supermultiplet. As was already mentioned, this supermultiplet includes the 1/2 BPS 
scalar operator O 20 ' j all the conserved currents and the Lagrangian of the theory. These operators 
appear as various components in the expansion of the supercurrent T in powers of the Grassmann 
variables. 


2.1 Stress-energy supermultiplet 

For 0 = 0 = 0, the lowest component of T is the scalar operator O^o'- It has conformal weight 2, 
belongs to the representation 20 ' of the 12—symmetry group SU (4) and has the following form 


O 20 '{x,y) = 0^^’^^ix)YABYcD = T{x,9,9,y) 


( 2 . 1 ) 


where the SU{4:) indices take values, e.g., A = 1,...,4. The auxiliary tensors YabYcd have 
been introduced to project the operator 0^^’^^{x) onto the representation 20'. They satisfy the 
relations Yab = —Yba and e^^'^^YabYcd = 0 and can be parameterized as 


YaB — ^cd^B — 


^ab Uab' 
Vba' y^f^a'b' 


( 2 . 2 ) 


where = det||yaa' ||= and we used composite indices A = (a, a') with a, a' = 1, 2 

and similarly for B = {b,b'). The variables yaa' have the meaning of coordinates of the operator 
O 20 'ix,y) on the four-dimensional coset SU{Y)/{SU{2) x 51/(2)' x 17(1)) of the 12—symmetry 
group. 

In addition, the supercurrent T describes a short supermultiplet of the 1/2 BPS type. This 
means that it is annihilated by half of the super-Poincare generators and, therefore, it effectively 
depends on 4 chiral and 4 antichiral Grassmann variables, 0“ and 0“, (with a, d = 1, 2 and a, a' = 
1,2), respectively. The stress-energy tensor can be obtained by applying A/ = 4 supersymmetry 
transformations to (2.1) and can be extracted from T with the help of the differential operator [11] 


T 


(x) 


{d9Ud9)ya{deha'{der- - 

+ ^{dx)iaddx)fj)pidy)aa'{dyT''" T(x, 0, 0, y) 


9 = 9=0 


(2.3) 


where we switched to spinor notation, (see Appendix A for details) and 

denoted the weighted symmetrization of the indices by, e.g., = \{Aap+A0a)- The hrst term 
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in the brackets in (2.3) selects the {9‘^)ap{d'^)ap component of the supercurrent, whereas the second 
and third terms involve total derivatives acting on lower components of T. These are typical 
conformal and symmetry descendant terms, which have to be subtracted in order to ensure 
the expected properties of the stress-energy tensor: current conservation = 0 

and zero i?—charge or, equivalently, independence on the y—variables, {dy)°' — 0- 

In a similar manner, the i?—symmetry current is related to the supercurrent T by the fol¬ 
lowing differential operator [11] 


Jad,aa'i^}y) (^5)00'(^6)oa T ^{dx^ad^dy^acir 


T{x,e,0,y) 


e=e=o 


(2.4) 


where the second term in the brackets is again the subtraction of a descendant needed to ensure 
the current conservation, Jad,aa'{x,y) = 0. 


2.2 Auxiliary spinor variables 

Let us start with the two- and three-point correlation functions of the supercurrent T(i) = 
T[xi,9i,6i,yi). The AA = 4 superconformal symmetry fixes them up to a normalization constant 
[23-27] 

(r(i)r(2)) = |(z)l2)^ 

(r(i)r( 2 )r( 3 )) = cDi2D2zD^i , (2.5) 


where Dij is a free (super) propagator 


y- y- 
n — ^ ^ 

- -2 “ ’ 

-^b 


( 2 . 6 ) 


with = xf^ -I- 


aa 

Hj 




fiOLa 

OLOL^ij 


for Xij = Xi — Xj and similarly 


iyi/)a'a9f^ and yfj^' = y\ 
for yij, 9ij and 9ij (see Appendix A for our conventions for rising and lowering Lorentz and the 
SU{2) indices). In the Af = 4 SYM theory with the SU{N) gauge group we have c = — 1. 

Putting 9 = 9 = 0 on both sides of (2.5) we reproduce the known expressions for the two- 
and three-point correlation functions of the 1/2 BPS operators (2.1). To obtain the correlation 
functions of the stress-energy tensor from (2.5), we have to apply the differential operator defined 
in (2.3) to each T(f) in (2.5). After a lengthy calculation we arrive at the expected result for the 
two-point correlation function (with xi 2 = xi — X 2 ) 


irr I \rr i \\ (^ 12 )(Q:ia 2 (®12)/3i)/32 (^ 12 )(a 2 «i (®12)/32)/3i 

(^aiai,/3l/3l(^l)^a2a2./32/32(*2)) = 160 C-- 


(2.7) 


Likewise, we can use the second relation in (2.5) to reproduce the known result for the three-point 
correlation function of stress-energy tensors in an AA = 4 superconformal theory [5] . 

We can simplify (2.7) by contracting the Lorentz indices of the stress-energy tensors with 
auxiliary light-like vectors re/ (with re? = 0) 


T(i) = nfnrre^x.) = jAfAf AfAf , 


( 2 . 8 ) 
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where nf“ = = AfAf with Aj and A* being arbitrary two-component (anti)chiral com¬ 

muting spinors. Then, we find from (2.7) 


(r(l)r(2)) = 160c 


((1|xi2|2](2|x2i|1])^ 


(2.9) 


where we used the standard notation for {i\xij\j] = Xf{xij)aa^'j ■ 

The auxiliary variables A^ and A* serve two main purposes. They automatically symmetrize 
the correlation function with respect to the chiral and antichiral Lorentz indices and, in addition, 
they simplify the conformal properties of the stress-energy tensor. To see this, we exploit the 
freedom in defining A* and A* to assign to them weights under conformal inversion I[{xi)adi] = 
{x^^)aa- Then, we choose [28] 


I[Xf] = w{xi){x-^ff^K0 , I[Af] = w{x,)\^{xry- , (2.10) 

with a scalar weight factor w{xi) = (x?)"" that can be chosen at our convenience. In this section, 
we put w{xi) = 1 for simplicity but, as shown in section 4.1, the choice of w{xi) = xj is more 
convenient for discussing the conformal properties of the four-point correlation functions. 

We would like to mention that the variables A* and A* have a natural interpretation in the 
context of scattering amplitudes where they are used to describe the helicity states of scattered 
massless particles. The invariance of the scattering amplitudes in planar A/” = 4 SYM under 
transformations (2.10) led to the discovery of the dual conformal symmetry in this theory [28]. 
As we will see in the next section, the final expression for the four-point correlation function 
of stress-energy tensors involves the same A/” = 4 dual conformal invariants that have been 
encountered in the study of the scattering amplitudes. 

With the assignment (2.10) and with w{xi) = 1 the quantity {i\xij \j] entering (2.9) transforms 
covariantly under inversion, /[(fjxjjjj]] = {i\xij\j]/{x‘^x‘j) leading to the simple transformation 
property of the two-point correlation function (2.9) 


/[(r(l)T(2))] = {xlxlf{T{l)T{2)) . (2.11) 


We deduce from this relation that the projected stress-energy tensor (2.8) transforms under 
conformal transformations as a scalar conformal primary operator with conformal weight 2, 
I[T{i)] = {xf)‘^T{i) [29, 30]. Notice that, by definition, T{i) is a homogenous polynomial in 
Xi and Xi of degree 2. These two facts will play an important role in what follows. 


2.3 Four-point correlation fnnctions 

In distinction with (2.5), the four-point correlation function of the supercurrents is not fixed by 
AA = 4 superconformal symmetry. The corresponding Ward identities imply that this correlation 
function can be decomposed into a sum of two terms carrying different information about the 
states propagating in various OPE channels 

Gi = (r(i)r(2)r(3)r(4)) = ( 2 . 12 ) 

Here receives contributions only from the protected operators in the theory whereas all 
unprotected operators contribute to , As a consequence, the form is fixed by A/" = 4 
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super conformal symmetry and is given by an expression analogous to (2.5) 


— c{Di2D23D34D4i + D 13 D 23 -D 24 -D 14 + .C)i2L>24.C)34L)i3) 

+ — (Z)^2-^34 +-^13-^24 +-^ 14 -^ 23 )) (2.13) 

where the second line corresponds to the disconnected contribution. 

The second term j^as more complicated form and, in contrast with (2.13), it is not a 

rational function of the distances x'fj. It admits the following representation 

g(anom) ^ (^Di 3 D 2 Af Mx, 6 , 6 , y) , (2.14) 


where the propagators Dij were defined in ( 2 . 6 ) and X 4 depends on the four superspace points. 
The product of propagators (Ili 3 T) 24 )^ carries the (super)conformal weights of the supercurrents 
whereas Z 4 is invariant under A/" = 4 superconformal transformations, X 4 = 0 with the genera¬ 
tors J = {Qo, 5^, Qi, S^,...} given in (2.20) below. 

As was shown in [11], the invariant X 4 admits two equivalent representations 


X4 = 




zj4 /j4 /j4 n4: 
^2 ^3 ^4 

^4 nA nA nA 
“1 “2 “3 ^4 


F{x) 

{Di^D2a)\ 

F{x) 

{Di^D2iY\ ’ 


(2.15) 


where = UaAQa and 5® = Od, denote the products of the 8 chiral Poincare super¬ 
charges and the 8 generators of antichiral special superconformal transformations, respectively, 
and similarly for the generators of opposite chirality. Since the generators Q and S are nilpotent 
and form an abelian subalgebra, {Q,<S} = 0 , it immediately follows from the first relation in 
(2.15) that Z 4 is annihilated by these generators. Similarly, it follows from the second relation in 
(2.15) that Z 4 is also annihilated by Q and S. 

The expression in the brackets in (2.15) contains the product of all the available (anti) chiral 
Grassmann variables, 9f = ][][^ ^ and Of = ][][q a' well as a scalar function F depending 

only on the coordinates Xi. The Bose symmetry of (2.14) implies that F{x) should be invariant 
under the exchange of any pair of points Xi Xj. In addition, the invariance of (2.15) under 
conformal transformations leads to the following inversion property of F [11] 

/ [F(x)] = ixlx2xixl)^F{x). (2.16) 


It allows us to express F(x) in terms of a function of conformal cross-ratios 

^{u, v) 


F{x) = 


u = 


■^12-^34 

™2 „2 ’ 
X 13 X 24 


V = 


2 

■■'^23-^14 

^2 2 
X 13 X 24 


(2.17) 


uv{xl^xl^Y ’ 

The symmetry of F[x) under exchange of points translates into the crossing symmetry relations 


<I>(u, v) 


$(u,u) = 

V 




(2.18) 


Substituting (2.15) and (2.17) into (2.14) we find that the general expression for 
depends on a single function ^(u,v). The same applies to the variety of four-point correlation 
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functions that appear as components in the expansion of powers of the Grassmann 

variables. In particular, putting = 0 on both sides of (2.14) we find that the fnnction 

<h(u,n) defines the four-point correlation function of 1/2 BPS operators O 20 '- The latter corre¬ 
lation fnnction has been thoronghly stndied in A/" = 4 SYM [31, 32] and the fnnction $(n, n) is 
known in this theory both at weak [33-36] and at strong conpling [37-39]. 


2.4 Master formula 


Let ns simplify relation (2.14). To begin with, we rewrite it using the first relation in (2.15) 


^(anom) 

^4 


(T>13T>24)" 


a4 n4 
^1 ^2 ^3 "4 


F{x) 

{DnD24? 


(2.19) 


where we used the composite index A = (a, a') to split the product of 8 generators into 
= Ha, a Qa and = Ha.' Qq and similarly for the product of generators and . These 
generators act on the snperspace coordinates (xj, 6i, Oi,yi) (with f = 1,..., 4) at each point and 
are given by the sum over the 4 points of the following anticommuting differential operators 


Qa = 


d 

de?: ’ 


q . _ _a 9 ^ 


~ ^ay^b'a o + Db' 
UXrvty 


d 


d 


QOaJ^aP ^a'yVb'Qya+^b'-^a'y 


d 


dB„ 


( 2 . 20 ) 


The explicit form of the A/” = 4 transformations generated by these operators can be fonnd in 
Appendix B. 

We notice that the snper-propagator (2.6) is annihilated by all the generators (2.20) except 
5? . This allows ns to simplify (2.19) as 


^|anom) ^ g4g/4^4y4 ^^4 ^4 ^4 ^4 ^ 


( 2 . 21 ) 


where the notation was introdnced for 


-2 


= {D^^D24YS^,^{P^^D24) 

d 

+ y^b' 




d 


t,aa 




^i,ay ^i,a‘ 




d 


dyf. 




d 




( 2 . 22 ) 


The only difference as compared with the last relation in (2.20) is that, in the first term, 
appears to the right of the derivative. 

As follows from (2.20), the generators Q“ and do not involve spatial derivatives d/dx 
and, therefore, in the expression on the right-hand side of (2.21), they only act on the prodnct of 
Grassmann variables. In this way, we arrive at the following relation ^ 


^(anom) 

^4 


[(a:i2a;i3^i4)^(a;i3^6'i3 - x^^ei2Y{x^lei4, - x^^ei2fF{x)\ , 


(2.23) 


f 


^To evaluate Q^S^(0f ... 9\) it is convenient to use the integral representation Q'^Oi 
and similarly for and apply relations (B.2) and (B.4) from Appendix B. 


f d‘^eexp(Q-€)F = 


7 















where the expression inside the brackets is invariant under the exchange of any pair of points and 
is independent of yi and 9i. This is the master formula that we shall use in the following section. 

According to (2.23), the dependence of on yi and 6 i comes entirely from the product 

of generators Q'^S'^ = WaaQZ' ^b' Obviously, the expression on the right-hand side of 

(2.23) is not symmetric under the exchange of chiral and antichiral sectors which seems to be in 
a contradiction with the expected reality property of the supercurrent T [40]. Nevertheless, this 
symmetry is restored owing to relation (2.15). If we use the second relation in (2.15) we would 
derive another equivalent form of (2.23) in which the variables 6 i and 6 i are exchanged. 


3 Four-point correlation function of stress-energy tensors 


In the previous section, we presented the general expression for the correlation function of the 
supercurrents, Eqs. (2.12), (2.13) and (2.23). To obtain the four-point correlation function of the 
stress-energy tensors 

GttttW = (3.1) 

we have to go through two steps: expand the supercorrelator t /4 in powers of Grassmann variables, 
and then apply to Q 4 the four differential operators defined in (2.3), one for each stress-energy 
tensor. 

As follows from (2.12), the resulting expression for this correlation function can be decom¬ 
posed into a sum of two terms, 

Gtttt = G^^rprp(x) -|- G^^^^(x). (3-2) 

Here the first term comes from (2.13). It is a rational function of the distances xjj, inde¬ 

pendent of the variables y*. Like the three-point correlation function, it can be decomposed into 
the sum of three different Lorentz structures coinciding with the four-point correlation functions 
of stress-energy tensors in a free theory of scalars, fermions and gauge fields, respectively. To 
save space, we do not present their explicit expressions. 

Let us turn to the second term on the right-hand side of (3.2) generated by (2.23). Since the 
differential operator in (2.3) is given by a sum of three terms, the operator that we have to apply 
to in order to extract contains 3^ = 81 different terms. Nevertheless, as we show 

in this section, it is sufficient to examine only one term, the one containing the maximal number 
of derivatives with respect to the Grassmann variables. It comes from the hrst term inside the 
brackets in (2.3) and selects the following component of the supercorrelator 


^(anom) 

^4 




ai/ 3 i...04/34 
oi/3i... 04/34 


+ . . . , 


(3.3) 


where (0^)“^ = ~ and dots denote the remaining components. 

3.1 Maximal number of derivatives recipe 

Let us first compute Matching (3.3) to (2.23), we observe that the 8 chiral variables 

9i on the right-hand side of (3.3) can only come from the expansion of the expression inside the 
brackets in (2.23). Moreover, since the latter expression has the same degree 8 in the 0’s, we can 


neglect the terms containing derivatives in the product of generators Q'^S'^ in (2.23). In this 
way, from (2.20) and ( 2 . 22 ) we obtain a simplified form of the generators Q' and S', 


d 


Qa' = Y.ha^ 

i 


^ dx • ^ 

d 


“ 1 “ „! 


d 




(3.4) 


Here in the second relation we also take into account that {—0^'^ Hb' does not contribute due 
to the y—independence of the expression inside the brackets in (2.23). 

In order to arrive at (3.3), we have to extract the (<9f)^^^^ ... component from the 

product of generators in (2.23). This turns out to be an extremely nontrivial task, mainly 

due to the presence of the second (nonlinear) term in the second line of (3.4). There is however 
a shortcut that allows us to overcome this difficulty. In what follows we shall refer to it as the 
‘maximal number of derivatives’ recipe. 

Let us examine the action of on a test function f{x, 6 ) independent of the variables 

Oi- Replacing the generators by their explicit expressions (3.4) we find that the expansion of 
Q"^S"^f{x, 6 ) involves various terms including those containing the maximal number (equal to 
8 ) of derivatives with respect to Xi. To identify such terms we can safely neglect the second 
(nonlinear) term in 5^,^ and, in addition, ignore the noncommutativity of Xi and .. Introducing 
the notation pi = dx^ and treating Xi and pi as commuting variables we get 


Q'^S^ 


^ 4 

f{x,0) = (X] 


2=1 

1 

& 


2=1 


E 


f{x,e) +... 

[^ 0 iPiXi 4 ^ f{x,9) + ... , 


(3.5) 


where we used the shorthand notation (OiPi)'^, = 9i^a'aPf^ and iOiPiXi)^,^ = Oi^ 0 i^pf°^x^ and 
denoted by dots terms with fewer spatial derivatives. By construction, (3.5) only describes the 
terms containing 8 spatial derivatives acting on a test function. 

We are now ready to formulate the maximal number of derivatives recipe. It consists of 
two steps: (i) expand (3.5) in powers of 9i and move all pi to the left of Xi, and then (ii) replace 
Pi = dxi so that the derivatives act on all Xi to the right. It turns out that the resulting expression 
gives the exact result for f {x, 9), including the terms shown by dots in (3.5). The proof of 

this statement can be found in Appendix B. 

Let us now apply the above recipe to identify the • • • (^ 4 ) 0 , 4^34 component of Q'‘^S'^f{x, 9). 

We hrst use the second relation in (3.5) to get 

Q'^5^/(x,6») = + ... , (3.6) 

N=i ^ 

where Maii 3 i...a 4 i 34 {x) is symmetric under the exchange of the spinor indices a* -H- /Ij. It admits 
a concise representation in terms of the auxiliary spinors introduced in ( 2 . 8 ) 

> ( 3 . 7 ) 
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where A4{x, A) is given by 


M{x,X) = (^(I|a;^2^x24|4)(3|a;3ixi2|2) - (I|xi2^x23|3)(4|x4ia;i2|2)^ , (3.8) 

with (l|x^2^a:24|4) = A“(x^ 2 ^a: 24 )o^A 4 ^/ 3 . Obviously, the expression on the right-hand side of (3.7) 
does not depend on the auxiliary spinors A,. They were introduced to make the symmetry of 
(3.7) under the exchange of indices more transparent. It is straightforward to verify that A4{x, A) 
is invariant under the exchange of any pair of points, (xj,Ai) -f-)- {xj,Xj). In addition, applying 
(2.10) with w{xi) = 1 we find that it transforms covariantly under inversion 

/[A4(x,A)] = {xlx2xlx‘l)~‘^M(x, X). (3.9) 


Let us now go through the second step of the recipe. Since all pi in (3.6) are located to the 
left of some x—dependent functions, we replace pi = dxi on the right-hand side of (3.6) and obtain 
the exact expression for this particular component of Q'^S^f{x,6). Notice that in the resulting 
expression the spatial derivatives act both on M.aii 3 i...a 4 i 34 ix) and on the test function f{x,9). 

Going back to (2.23) we choose the test function in (3.6) to match the expression inside the 
brackets in (2.23), 

f{x,e) = {xl2xl3xl^f{x];^ei3 - X^ 2 ^l 2 )^{x^leu - X^2^12)'^F{x) . (3.10) 


In order to compute (3.3) we only need its component proportional to (0f)ai/3i ■ • • (^ 4 )a 4 / 34 - Going 
through the calculation we find 


f{x,e) = ]J{efr^^M^,f,,...^,f,,{x)F{x) + ... , (3.11) 

i 

where ■M.aii 3 i...a 4 ,p 4 ,{x) is given by (3.7) and the dots denote other components. Finally, we 
substitute the last relation into (3.6) and match the result to (3.3) to find 


G 


ai/3i...04/34 
04/34 


i=l 


•^^oi/3i... 04/34 (^).Ad.Yi(5i...74(54 i.X^Fi^x') 


(3.12) 


Here the nontrivial information about the particular A7 = 4 superconformal theory is encoded in 
the function F{x) defined in (2.17). 


3.2 Restoration of the symmetry 

We recall that relation (3.12) defines a particular component of the supercorrelator (3.3). To get 
the four-point correlation function of the stress-energy tensors (3.1), we have to add to (3.12) the 
contribution of the remaining 80 terms mentioned in the beginning of this section. To understand 
their role, let us check whether (3.12) is consistent with the expected properties of the stress- 
energy tensor. 

To simplify the analysis, consider the correlation function 

= ( 3 . 13 ) 
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where the dots denote other operators. Then, the symmetry and conservation properties of the 
stress-energy tensor, = 0, yield the Ward identities ^ 

_ ai/3i ^oi/3i _ /o \Qi^«i^i _ r. /r. .. .x 


In addition, should transform covariantly under conformal transformations with weight at 

point xi corresponding to a conformal primary operator of Lorentz spin (1,1) and dimension 4. 

Let us verify whether (3.12) satisfies relations (3.14). It is easy to see that the first two 
relations in (3.14) are automatically satisfied due to the symmetry of (3.7) under the exchange of 
indices, M.aii 3 i...a 4 i 34 {x) = 0. Examining the last relation in (3.14) we find, using the identity 
dada'y = 


(ri i^^l$l---a4$4 

04/34 


gai7lQ 


XI 


•^^ai/3i...ci4/34 (^)-^^7 i5i...74(54 (^)'^(^) 


(3.15) 


Since this expression is different from zero for a generic function F(x), we conclude that (3.12) 
does not respect the conservation of the stress-energy tensor. This is not surprising since (3.12) 
is only a part of the correlation function (3.1) and for the expression on the right-hand side of 
(3.15) to vanish we have to add the contributions of the remaining 80 terms. 

We notice that (3.15) would automatically vanish if the expression inside the brackets on 
the right-hand side of (3.15) and (3.12) are symmetric under the exchange of indices ai and 71 . 
Similarly, the vanishing of would follow from the symmetry under -H- 7 ^. 

This suggests that the net effect of the remaining 80 terms is to symmetrize the product of two 
tensors in (3.12) with respect to the indices (a*,/3j, 7 *, <5*). Taking into account (3.7) we find 
that (up to an overall normalization factor) this amounts to replacing in (3.12) 


A^«i/3i...a4/34(a^)-^7i5i...74<S4(^) —^ (^A^) 7 i [7W(x, A)]^ (3.16) 

i 


with A4(x,A) given by (3.8). Upon this substitution, (3.12) respects the conservation of the 
stress-energy tensor. 

Let us now examine the conformal properties of (3.12). As before, it is convenient to examine 
the inversion. Denoting the expression inside the brackets in (3.12) by Fa/s^s, we find using (2.16), 
(3.7) and (3.9) that it transforms covariantly 






(3.17) 


Then, we use = 0^=1 dxVdxf"Fap^^ and apply the identity (A.5) to find 


Tr^Ol/ 3 i...04/34 
1^Q1;3i...04/34 



(3.18) 


^Strictly speaking, relation (3.14) is valid up to contact terms [2, 41, 42]. For four-point correlation 
functions in an A/” = 4 superconformal theory, the contact terms are only due to the rational part of (3.2), 
whereas the anomalous contribution to (3.2) is less singular and satisfies homogenous Ward identities. 
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where the dots in the second relation denote inhomogenous ‘bad’ terms coming from the commu- 

1 -J.„ ---- j.—^ai/3l...d4/34 

04/34 


tators [dxidx\ \ ]. Due to the presence of such terms, does not transform 


covariantly. Their contribution however involves terms like e" and it would vanish if 

F„//3'7 IS' is symmetric with respect to the indices a' and 7 '. This is exactly what happens if we 
apply (3.16) to (3.12). 

Thus, we conclude that the substitution (3.16) not only ensures the conservation of the 
stress-energy tensor but also restores the correct conformal properties of the four-point correlation 
function. 

The question remains however whether (3.16) correctly describes the additional contribution 
to (3.12) coming from the last two terms on the right-hand side of (2.3). Going through a lengthy 
and tedious calculation we verified that this is indeed the case. 


3.3 Simplified form of the correlation function 

Combining together (3.12) and (3.16), we arrive at the following remarkably simple result for the 
four-point correlation function of stress-energy tensors 

4 

1 = 1 


More precisely, this relation describes the ‘anomalous’ (non-rational) contribution to (3.2). We 
recall that F{x) is the nontrivial dynamical function defined in (2.17) and A4{x, A) is a kinematical 
function given by (3.8). Relation (3.19) is one of the main results of this paper. 

Since A4{x,X) is a homogenous polynomial of degree 2 in A", the expression on the right- 
hand side (3.19) does not depend on the auxiliary spinors Aj. This polynomial admits an elegant 
representation when expressed in terms of new variables 


= 


(xO^^A, 


i/3 



I = (ce, d). 


(3.20) 


which satisfy the following determinant relation 


{1234) = eijKLZizizi^Z^ 

= (l|x5'2^X23|3)(4|x4ia;i2|2) - (I|x5'2^x24|4)(3|x3ixi2|2) . (3.21) 


Then, we find from (3.8) 


M{x,\) = (1234)2. 


(3.22) 


In this representation, the properties of A4{x,\) becomes more transparent. 

Relations (3.20) and (3.21) are very familiar from the dual space description of scattering 
amplitudes. There, the variables A* and pi = Xi — Xj+i define the helicity and the on-shell mo¬ 
mentum of each scattered particle, respectively, and the dual conformal symmetry is realized 
as SL{4) transformations of Z-. Relation (3.21) defines the simplest four-point dual confor¬ 
mal invariant and serves as a building block in constructing four-dimensional integrands for the 
scattering amplitudes in planar J\f = 4 SYM [43, 44]. 

As explained in the previous subsection, relation (3.19) respects the conservation and confor¬ 
mal symmetry of the stress-energy tensor, independently of the form of F{x). The expression on 
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the right-hand side of (3.19) involves 8 spatial derivatives acting on the product of two functions. 
Obviously, its expansion yields a very lengthy expression involving the 22 different Lorentz struc¬ 
tures mentioned in the Introduction, each of them involving derivatives of F{x). The very fact 
that the correlation function in an A/” = 4 superconformal theory admits a compact representation 
(3.19) leads to important consequences that we explain in section 5. 

We can further simplify (3.19) by considering the correlation function of the operators (2.8). 
This amounts to projecting the Lorentz indices on both sides of (3.19) with the auxiliary spinors 
Aj and A* 


G 


(anom) 

TTTT 




04/34 

'-^oi/3i... 04/34 


i=l 


(3.23) 


Introducing the notation for the differential operator 


= Xi,a{9xiT°‘{dxi)a = [i\dxi\dxi) 


(3.24) 


we finally obtain from (3.19) a very compact representation for the correlation function 

= 4^ (T>iT> 2 T> 3 T' 4 )^ [(1 2 3 4)^F(x)j. (3.25) 


Here the additional normalization factor comes from = 4(12 34)^. 

We notice that, at first glance, relations (3.19) and (3.25) are not symmetric under the 
exchange of the chiral and antichiral sectors, a* •H'd* and Pi ^ pi. As was already mentioned in 
the previous section, thanks to relation (2.15), the correlation functions (3.19) and (3.25) admit 
another equivalent representation in which the chiral and antichiral indices are interchanged. 


4 Four-point correlation functions of currents 

In the previous section we demonstrated that the calculation of the correlation function of stress- 
energy tensors can be greatly simplified by employing the maximal number of derivatives recipe. 
It allowed us to obtain the very special representation (3.12) for the relevant component of 
the supercorrelator ^ 4 , and then to promote it to the complete expression for the correlation 
function by simply symmetrizing the product of Lorentz tensors with respect to the chiral indices, 
Eq. (3.16). 

In this section we argue that the same approach can be applied to the four-point correla¬ 
tion functions involving other components of the stress-energy supermultiplet, for example, the 
i?—symmetry current (2.4). 

4.1 Ansatz 

Instead of going through the computation of the various correlation functions, let us try to 
generalize (3.19) and formulate an ansatz for the correlation function involving conserved currents 
of Lorentz spin (5/2,S'/2) 

GS:::“f(x,...) = (o|j“\-:4nx)...|o), (4.i) 

where the dots denote the remaining operators. For 5 = 1 and 5 = 2 the operator (^) 

coincides with the ii—symmetry current and the stress-energy tensor, respectively. In close 
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analogy with (3.12), we assume that the correlation function (4.1) has the following general 
form [11] 



g«171 _ _ _ QOis'rs 


M 


«171---“S7S 


[x)F{x) 


(4.2) 


where dx'^ = djdx-ya and the number of spatial derivatives matches the Lorentz spin of the 
operator. If other operators in the correlation function (4.1) carry nonvanishing Lorentz spin, 
the expression on the right-hand side of (4.2) contains additional derivatives acting on their 
coordinates. 

To ensure the conservation of the current, we require that M.ai'Yi.-.as'rsi ^^...) be completely 
symmetric with respect to all 2S chiral indices. Indeed, we verify that 






M 


Ql7i...0:575 


ix)F{x) 


= 0 


(4.3) 


for an arbitrary x—dependent completely symmetric tensor ■Mai'yi...as'ysi^^ ■ • 
demand that (4.2) has the correct transformation properties under inversion ^ 


■ ■ ■ )l = ■ ■ ■), (4.4) 

where we do not display the weight factors corresponding to the other operators. Taking into 
account (2.16) and the identity (A.6), we find that (4.4) and (4.2) lead to 

/[A^ai7i...as7s(a:)] = Mfj,Si...fSs 6 s{x) ■ (4.5) 


As before, we can simplify this relation by projecting all chiral Lorentz indices by the auxiliary 
spinor A, 


Ms(x,A) = (4-6) 

Taking into account the transformation properties of the auxiliary spinors under inversion (2.10), 
we find from (4.5) 

I[Msix, A)] = w‘^^{x){x‘^)~^~‘^Msix, A), (4.7) 

where w{x) = (x^)” with an arbitrary n. In the previous section we used w = 1 but it is now 
more convenient to choose w = x^, so that /[A"] = x"^A^. 

The X— dependent factor on the right-hand side of the last relation defines the local conformal 
weight associated with the current of spin S. If the correlation function (4.1) contains four currents 
with different spins Si, the total conformal weight of the corresponding function AI 5152 S 3 S 4 (^) is 
given by the product of four such factors, one for each current. As a consequence, A45j5253S4(®) 
depends on four auxiliary spinors Aj (with i = 1,... ,4) and satisfies two main requirements: (i) 
to be a homogenous polynomial in Xf of degree 2Si and (ii) to transform under inversion as 

^[-ii^Si52S3S4(®) i^)] = AIsi 525354 ( 3^5 ■ (4-8) 

i 

In addition, if two operators are identical with, e.g., 5i = 82 , in virtue of the Bose symmetry of 
the correlation function, the function A45^5^5354(x, A) should be invariant under the exchange of 
points 1 and 2. 

"‘We recall that inversion swaps the chirality of the Lorentz indices. 
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Having determined A^S’iSi53S4(ic, A), we can reconstrnct the fonr-point correlation fnnction 
of the cnrrents Jal...as^ (xi). To simplify the resnlting expression, we project all Lorentz indices 
by the (anti)chiral auxiliary spinors Aj and Aj and define 


^<554 T 


Ga&aa(x,A,A) = A? ...A7ni4. (4.9) 

This function admits the following representation (up to an overall normalisation factor) 


G5iS253S4(a^,A,A) 


(4.10) 


where the differential operator Vi is defined in (3.24). This relation generalizes (3.25) to the 
case of currents of spin Si = 1,2. Moreover, it is also applicable for S'* = 0 in which case the 
corresponding operator is the 1/2 BPS operator O20'- Notice that the i?—symmetry current 
and the 1/2 BPS operator have an i?—charge and depend on the auxiliary variables y. In the 
representation (4.10), this dependence is carried by the function A45^525354(x, A). 

As shown in [11], relations (4.2) and (4.10) can be generalized to the fermionic A/" = 4 super- 
symmetry currents Jaia 2 i which are also members of the stress-energy supermultiplet. 

They carry Lorentz spins (3/2,1/2) and (1/2, 3/2), respectively, and should appear in pairs in a 
non-vanishing correlation function. The only difference as compared with (4.2) is that the numbers 
of (chiral) a— and 7—indices of the fully symmetric tensor M do not match. As a result, in the 

general case of operators of Lorentz spin {Si/2, S[/2), their four-point correlation function is given 
S' s'' 

by (4.10) with V^ replaced by V^'. As already mentioned, there exists another, equivalent rep- 
resentation of the same correlation function involving the conjugate operators = [d\.\dxi\i)^' 
and conjugate antichiral tensor A45'5'5^5'(x. A). 


4.2 Special solutions 

Let us construct some solutions to (4.8). We notice that (4.8) defines At5^52S'3S’4(^) •^) to 
multiplication by an arbitrary function of the conformal cross-ratios u and v defined in (2.17). 
Such a function will in general induce additional singularities for xfj —>-0. To fix the ambiguity, 
we shall assume that A453525354(^1 A) should not have such singularities. 

We expect that A45^52S3S4(^) A) is a rational function of the distances admitting an 
analytic continuation to complex space-time coordinates. The rational behind this is that the 
analysis of the conformal properties of Alsj52S'3S4(^) •^) t)e simplified by employing Dirac’s 
embedding formalism. There the complexified Minkowski space is realized as a light-cone in 
complex projective space with homogenous coordinates = —X'^^ (with I,J = 1,... ,4) 
satisfying ejjklX^'^X^^ = 0. The complex coordinates dehne a particular parameterization 
of X^-^ 


xU 




Xij = -eijKiX^^ 


X ea/3 X^^ 


(4.11) 


with composite indices I = (a,Q;) and J = {13,(5). Then, the conformal transformations corre¬ 
spond to global SL{4;C) transformations of X^'^. The attentive reader will notice the similarity 
between (4.11) and (2.2). Indeed, Xaa and y^a' appear on an equal footing as bosonic coordinates 
of the supercurrent. 
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In addition to (4.11), we also need the variables (3.20) that carry the dependence on the 
chiral spinors A. As was already mentioned, there exists some freedom in choosing the weight 
factor w in (2.10). The advantage of the choice w = \s that the action of inversion on 
corresponds to a global 5'L(4) transformation ^ 

= 0SZ/. (4.12) 

We can use the variables and Zf (with i = 1,...,4) to define various 5'L(4) invari¬ 
ant quantities. Taking into account the identity Z:[{Xi)jj = 0 we can define three different 
Z—dependent structures 

{1234:) = eijKLZiZ^Z^Zi:, 

'^[12]3 = Z{{Xz)ijZ2 ~ (l|xi3a;32|2), 

Ai[234] = Zi{X2)lj{X^Y^{X^)KL{Zl)^ ~ (1|xi2X23X34X4i|1) , (4.13) 

where the first structure already appeared in (3.21). X^i2]3 A’i[ 234 ] are antisymmetric with 

respect to the points indicated inside the brackets, whereas (12 3 4) is completely antisymmetric 
with respect to the four points. We verify that, as expected, they transform covariantly under 
inversion 

I[(12 34)] = (12 3 4) , I[A’[i 2 ] 3 ] = — A’[ 12]3 5 2 2 2'^1[234]- (4-14) 

X 3 X2X'^X^ 

We are now ready to construct solutions to (4.8). 

Let us first revisit the four-point correlation function of stress-energy tensors, 5i = S '2 = 
Sz = Si = 2. Denoting the corresponding solution to (4.8) as A4tttt{x, X), we deduce from 
(4.8) that it should take the form of a homogenous polynomial in Zf^ of degree 4, invariant under 
inversion and under the exchange of any pair of points. Examining (4.14) we immediately find 
the solution 

X\tttt{x^X) ~ (12 3 4)^, (4-15) 

which agrees with (3.25). There exist other solutions, e.g. <Ti[234]‘T2[34i]A:’3[4i2]<T4q23]/ni<i 
but in contrast with (4.15) they have singularities for a;?- 0. 

The second example is the correlation function of two stress-energy tensors and two 1/2 
BPS operators, Si = S '2 = 2 and S 3 = S 4 = 0. The corresponding function, A4ttoo(x, X), is a 
homogenous polynomial in Zi and Z2 of degree 2 . According to (4.8), it should have conformal 
weight l/(x|x|)^ under inversion and be invariant under the exchange of points 1 ■<->• 2 and 3 -H- 4. 
We use (4.14) to obtain 

Mttoo{x,X) ~ (A’[i2]3<T[i2]4)^ = [(IIX 13 X 32 12)(l|xi4a;4212)]^ . (4.16) 

As in the previous case, there are other solutions involving A’i[ 234 ] but we have to discard them 
since they are singular for xf^ —>■ 0. In distinction with (4.15), the function Mttoo should 
also depend on the variables y at points 3 and 4. This dependence is unambiguously fixed 

^In fact, can be identified as coordinates on twistor space. 
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by the i?—symmetry and amounts to an additional factor of (j/34)^ on the right-hand side of 
(4.16). This factor has the meaning of the Clebsh-Gordon coefficient corresponding to the singlet 
representation in the tensor product decomposition 20 ' x 20^ 

The third example is the correlation function of two stress-energy tensors and two symmetry 
currents, 5i = ^2 = 2 and 53 = ^4 = 1. Going along the same lines as in the two previous cases, 
we arrive at 


MTTJjix,X) ~ (12 34)2T[i2]3^[12]4 = (1 2 3 4)^ (1 1x13X3212) (1 |xi4X4212). (4.17) 

The y—dependence of the i?—current y*) leads to an additional factor of on the 

right-hand side of this relation. As in the previous case, it corresponds to the singlet representation 
in the tensor product 15 x 15. 

We notice that, up to an overall normalization factor, the obtained expressions (4.15), (4.16) 
and (4.17) satisfy the following interesting relation 

~ Mtttt -XAttoo ■ (4-18) 


4.3 Comparison with known results 

Combined with (4.10), relations (4.16) and (4.17) yield predictions for the correlation functions 
involving two stress-energy tensors. The same correlation functions can be computed using the 
approach described in section 3. 

As an example, we can consider the correlation function (T(1)T(2)0(3)0(4)) containing 
two 1/2 BPS scalar operators. Following (3.3), we identify the ni=i component 

of (2.23), and then apply the maximal number of derivatives recipe to obtain a representation 
similar to (3.12). In close analogy with (3.16), the restoration of the conservation of the stress- 
energy tensor and its conformal properties can be achieved by symmetrizing the corresponding 
A4—tensor with respect to the Lorentz indices. Going along these lines we obtain 


(T(l)r(2)0(3)0(4)) 


{V 1 D 2 Y (y|4)^(l|xi3X32|2)^(l|xi4X42|2)"'F(x) 


(4.19) 


In the same way, for the correlation function containing two 7?—symmetry currents we find 


(r(l)r(2)J“3“3(3)J“'‘“4(4)) 


= {ViV2Yv^V^ 


y3T"y3T"(1234)2(l|xi3X32|2)(l|xi4X42|2)F(x)l, 


(4.20) 


where J“*“i(f) = JYYYxi)XfXf. Gomparing these relations with (4.16) and (4.17) we observe 
perfect agreement. 

We can apply the same technique to computing the correlation functions of other components 
of the stress-energy supermultiplet. The components T = ... -|- 9'^L{x) + 9'^L{x) -|- ... are of 
particular interest since they dehne the chiral (L) and antichiral (Z) on-shell Lagrangians of the 
theory, with L[x) = L^x) and i{L{x) — L{x)) being a topological term. 


(r(i)r(2)L(3)Z(4)) = {v,V2Ynl^ [Mttll{x,x)f{x)] , 
(L(1)Z(2)L(3)Z(4)) = [{xl,fF{x)\ , 


(4.21) 






with Mrprp^j^[x, A) = = (l|a;i4a:42|2)^. It easy to verify using (A.7), (4.8) and (2.16) that 

the expressions on the right-hand side of (4.21) have the correct transformation properties under 
inversion. The second correlation function in (4.21) is equivalent to the result of [45] obtained by 
a different method. The latter has been used as a nontrivial consistency check of two AdSs x S® 
supergravity calculations, that of the dilaton/axion amplitude in [46] and that for the bottom 
component of the same massive multiplet in [38]. 

Another class of correlation functions containing three 1/2 BPS operators and a conserved 
current has been studied in [11]. We can use these results to obtain 

(r(l)0(2)0(3)0(4)) =p2 [yl^yl^yUX,y^^,^{x,\)fF{x)\ , 

(J-'(l)0(2)0(3)0(4)) = Pi[Ti%4](y>^)'^i[234](x, A)F(x)], (4.22) 

where Xi^ 234 ] was defined in (4.13) and the notation was introduced for 

= ®12^34y23y24^1[34] + 3;23®142/24y34^1[23] + 2:13x142/231/34^1 [42] j (4.23) 

with = {yiiyijyjiY°‘ ■ Here Ti[ 234 ] and carry the dependence on the //—variables and 
are completely antisymmetric under the exchange of the points indicated inside the brackets. We 
recall that the i?—current and the 1/2 BPS operator belong to the <S'C/(4) representations 15 and 
20', respectively. The three terms in the expression for Ti[ 234 ] correspond to three overlapping 
5C/(4) representations in the tensor products 15 x 20' and 20' x 20'. 

Thus, we demonstrated on various examples that the correlation functions of conserved cur¬ 
rents in an A/" = 4 superconformal theory have the general form (4.10) in which the information 
about the quantum numbers of the currents is encoded in the A4—function. In particular, if the 
currents carry i?—charge, the A4—function can be decomposed into a sum of various irreducible 
components corresponding to overlapping SU (4) representations in the different channels. 

We would like to emphasize that the expressions for the correlation functions presented in 
this section are valid up to overall normalization factors. The latter can be determined, e.g., from 
the consistency with the operator product expansion. 

5 Application to energy-energy correlation 

In this section, we use our results for the four-point correlation functions to compute the energy- 
energy correlation (EEC) in an AA = 4 superconformal theory. This infrared safe observable 
describes the flow of energy in the final state created from the vacuum by some source [12]. 
More precisely, EEC(ni,n 2 ) measures the correlation between the energy fluxes in two different 
directions defined by light-like four-vectors ni and n 2 (with nf = = 0). As such, it is expected 

to be a regular positive-definite function of the angles defining the relative orientation of ni and 

742. 

5.1 Generalized optical theorem 

The energy-energy correlation admits the following representation in terms of a correlation func¬ 
tion [13-17] 

EEC(ni, 77 . 2 ) = J G(x; ni, 712 ) , 

G{x 3 - X 4 ; 741 , 742 ) = {£{ni)£{n2)J{x3)J{x4))w , (5.1) 



18 


where the Fourier integral fixes the total momentum of the final state to be q. The normalization 
factor a is fixed by the requirement that, in the rest frame of the source, for q = (E, 0) and 
Hi = (1, nj), the energy-energy correlation averaged over the spatial orientation of the unit vectors 
Hi has to satisfy the condition f 712 ) = 1. 

The correlation function G(x 34 ;ni,n 2 ) involves the so-called energy flow operators, £(ni) 
and T(n 2 ), which are expressed in terms of integrated stress-energy tensors (see [14, 15] for the 
explicit expression). The operator J(x) dehnes the source. If the source is created by a conserved 
current of spin S, the operator J(x) depends on the polarization tensor e 


= (5-2) 

The operator J(x) = J^x) is given by a similar expression with e replaced by the conjugated 
polarization tensor e. The subscript w in the second line of (5.1) indicates that this is the 
Wightman (not time-ordered) four-point function. It is related to its Euclidean counterpart via 
analytic continuation. 

To compute (5.1) it proves convenient to perform a conformal transformation and go to new 
coordinates x^ —>■ [20, 47] 


,1 __ XX ^ X ,2 

z =-— , z = X -— , z = — , dz = 

x+ x+ x+ 

where the notation was introduced for the light-like coordinates x^ and x = 


(^’ 

(x,x), 


(5.3) 


Xaa — ^/x(*^^)od — 


X^ X 


X X 


dx^ = dx~^dx — dxdx . 


(5.4) 


Then, if the source (5.2) is defined by an operator of spin S and dimension A 5 , the function 
G{x 34 ]n,n') in (5.1) is given by the four-point correlation function integrated over the light-cone 
coordinates of the two stress-energy tensors 


Gsix34;ni,n2) = —TA- 


ai-o^s e^i-^s^«i7i/3w/3i5w4', A“®7S('3jA^s^S('4') 


(n+n+)3 di-Hs 

poo 

X / dz^dz^ (T(0+, zf, Ti)r(0+, 2;2', T 2 ) (^4))vF , (5.5) 


where the product of the 2 :—dependent factor and A“7(i) = dz'Jdxf^ in the first line arises from 
the conformal transformation (5.3) of the two currents. Here T{zf = Q,z^,Zi) (with i = 1,2) 
are the stress-energy tensors (2.8) whose Lorentz indices are projected onto the same auxiliary 
spinors Aq = (1,0) and Aq = (J) satisfying Aq^q^Aq = z^. The nontrivial dependence of (5.5) 
on the light-like vectors ni and n 2 comes through the variable Zi = Hijnf, depending on the 
light-cone coordinates of these vectors. 

According to (5.1) and (5.5), the definition of the energy-energy correlation depends on the 
choice of the source J. Later in this section, we consider four choices for J: a 1/2 BPS operator, 
an i?—symmetry current, a stress-energy tensor and a Lagrangian of the theory.® We will show 
that, in virtue of AA = 4 superconformal symmetry, the energy-energy correlations are given in 
all four cases by the same expression and, therefore, do not depend on the choice of the source. 

®In all cases except the last one we have As — S = 2. The Lagrangian appears as 0{9‘^) component of 
the stress-energy supercurrent T and it has dimension 4 and spin zero. 
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5.2 Integrated correlation functions 

Let us first examine the Euclidean version of the correlation function entering (5.5). In close 
analogy with (3.2), this correlation function can be split into a sum of rational and anomalous 
pieces. As was shown in [18, 19], the contribution of the former to the energy-energy correlation 
(5.1) vanishes for generic ni and n 2 and, therefore, can be discarded.^ 

The anomalous contribution to the correlation function in (5.5) can be found from (4.10) 
for 51 = S 2 and S 3 = S 4 = S. More precisely, in order to match (5.5), we have to identify the 
auxiliary spinors, Ai = A 2 = Aq and Ai = A 2 = Aq. Taking into account that Aq and Aq satisfy 
for an arbitrary four-vector z, we can simplify the expression for the differential 
operators and X >2 defined in (3.24) (see (A.4)) ® 


A = Ao,d(dzJ“"(9Aja 


-XQ{dx,)adsi 



(5.6) 


In addition, we can safely neglect terms containing d - = djdz- since they produce a contribution 
to the correlation function (4.10) that integrates to zero after substitution in (5.5). Since the 
expression inside the brackets in (4.10) is a homogenous function of Aj, we can replace T>i ~ 
for i = 1, 2. In this way, we obtain from (4.10) 


(r(i)r(2)j5(3)J5(4)) 


{dM\V3Vif 


MttJsJs{^^X)F{z) 


+ ... , 


(5.7) 


where the dots denote terms containing total derivatives with respect to z^ and Z 2 , and we used 
the shorthand notation for T(z) = T(0+,2;~,Tj) and 

Js{k) = XfXka, ... A^"Afc,.^,4i4|(zfc). (5.8) 


For a source defined by the Lagrangian of the theory, we hnd from (4.21) that the correlation 
function (5.7) takes a slightly different form 


(T(1)T(2)L(3)L(4)) 




■X^TTLhi^^ X)F'{z) 


+ 


(5.9) 


As before, to obtain the correlation function in (5.5), it suffices to differentiate both sides of (5.7) 
with respect to the auxiliary spinors Xk and Xk with A; = 3,4. 

At the next step, we use the expressions for MttJsJs obtained in the previous section to 
evaluate them in terms of the z—coordinates. For instance, for the 1/2 BPS operator O 20 ' we 
apply (4.16) to get (with Zij = Zi — Zj) 

Mttoo = (2/34)^[(01^132:32|0)(0|2;i 4242|0)]^ = (2/i4)^^i2(4■ (5-10) 


Here in the first relation we rewrote (4.16) in z—coordinates, added the y—dependent factor 
and identified the auxiliary spinors |1) = |2) = |0). In the second relation, we replaced {zi)aa 
by its expressions (5.3) and (5.4) in terms of light-cone variables and took into account that 
zt = z+ = 0. 

^More precisely, it is proportional to a delta function with support at x = {nin 2 )/{ 2 {qni){qn 2 )) = 1. 

^Notice that we have to apply to derivative d\. before putting Ai = Aq. 
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For the Lagrangian, the stress-energy tensor and the current, we apply relations (4.21), 
(4.15) and (4.17), respectively, and go through similar calculations to find 

^ttlL = ^ 12 (^ 4 ")^) 

Mtttt = ^i2(3|w|4)^ , 

Mttjj = ytT^y3T^^t24^4 (3|w|4)2 , (5.11) 

where (3|a;|4) = A 3 cua/s and the matrix cua/s is given by 


LO = 


0 4 

-Z4 Z34 


(5.12) 


Examining (5.10) and (5.11) we observe an interesting property: the obtained expressions 
for the Ad—functions are independent of the coordinates and and, at the same time, they 
have the same dependence on zu- The former property implies that, upon substitution of (5.7) 
and (5.9) into (5.5), the integration over z^ and can be reduced to evaluating the following 
integral 

1 

/ = .^2 ,2 > (5-13) 

J-00 1^12^341 ^3 ^4 

where the function F{z) was defined in (2.17). Here the 2 :—dependent factor on the right-hand 
side carries the scaling dimension of the integral, so that ^( 7 ) is dimensionless. Moreover, the 
detailed analysis shows [11, 18, 19] that the argument of ^( 7 ) is a rational function of 7 = ^{zi) 
whose form is much simpler when expressed in terms of x—coordinates 


_ ^ {X34ni){x34n2) 

xl4{nin2) 


(5.14) 


The explicit expressions for ^( 7 ) in AA = 4 SYM at weak and at strong coupling can be found in 
[18]. We do not need them for our purposes. 


5.3 EEC for a scalar source 


To avoid the technical difficulties related to the Lorentz structure of the source, let us first consider 
(5.5) with a scalar source J given by the 1/2 BPS operator O 20 ’ oi’ the Lagrangian L. We apply 
(5.7) - (5.11) and (5.13) to get for Ao^^, = 2 , = 4 and S = 0 


Go{x34-,ni,n2) 


GL(x34;ni,n2) 


iyl4? 






izi2f{d^^d^4[zi2Gh)], 


4)^0 


( 44 )" 


24 , 


(444244)^ 


( 212 ) [^;i 2 ^( 7 )], 


(5.15) 


where the subscript on the left-hand side indicates the choice of the source. We can now use (5.3) 
to write the expressions on the right-hand side in terms of Xj. This can be done with the help of 
the identity 


[ 7 ^( 1 -7)2^"(7)]" = {z42f{d^^dz4[zi2G{l)] = (5-16) 

ID X 34 
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which can be verified by a straightforward calculation. Taking into account that 2 (nin 2 ) = 
X 34 = 2 : 34 /we arrive at Lorentz covariant expressions for the integrated 
correlation functions 


Go(a:34;rai,n2) 


GL{xi,i]ni,n2) 


{yh? ^2 £( 7 ) 

128(nin2)3 xi4 ’ 

_ \ _^2 ^2 

128(nin2)3 xi4 ‘ 


(5.17) 


We notice that the two expressions are related to each, Gl Ul^Go. 

To compute the energy-energy correlation (5.1), we have to analytically continue relations 
(5.17) to get the Wightman correlation functions, and then Fourier transform them with respect 
to X 34 = X 3 — 2 : 4 . It is easy to see that the Fourier transforms of Go and Gl are proportional to 
the same dimensionless function, e.g. 


J d^x e“'^ Gl (x; ni, 77-2) 


128(nin2)3 


Hx), 


(5.18) 


where the notation was introduced for x = q‘^inin 2 )/{ 2 {qni){qn 2 )) and 


Hx) = q^ j 


x2 — zOx*^ 


(5.19) 


Here Qw{l) corresponds to the particular analytic continuation of its Euclidean counterpart (/(y). 
Together with the zOx*^’ prescription in the denominator this ensures that the Fourier integral 
(5.19) is different from zero only for > 0 and qo > 0, as it should be for a physical quantity 
measuring the flow of energy in a final state with total momentum q. As a consequence, J-{x) is 
different from zero for 0 < x < 1- In the rest frame of the source, for q = {E,0) and = (l,ni) 
(with nf = 1), the scaling variable x = (n.in 2)/2 = (1 — cos 6)/2 is related to the angle 9 between 
the unit vectors ni and n 2 - 

Substituting (5.17) into (5.1) we find that the energy-energy correlations in both cases are 
proportional to E{x)- The proportionality factor can be determined by imposing the normalisa¬ 
tion condition 

[\xEEC{x) = l, (5.20) 

Jo ^ 

which follows from the requirement for the total energy in the final state to be equal to the 
momentum transferred q. In this way, we arrive at 

EECo =EECl = (5.21) 

with X = (^^ 1 ^ 2)72 = (1 — cos0)/2 in the rest frame of the source. 


5.4 EEC for a tensor source 

For a scalar source, Lorentz invariance implies that EEC(ni,n 2 ) can only depend on the relative 
angle x between ni and n 2 - For a source dehned by the i?—current and the stress-energy tensor, 
this is not necessarily the case due to the dependence of EEC(ni,n 2 ) on the polarization vectors 
(5.2). Nevertheless, as we show in this subsection, the energy-energy correlations in an AA = 4 
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super conformal theory do not depend on the choice of the source and are given by the universal 
scaling function (5.21). 

For a source given by a current of spin S, an additional complication arises in (5.5) due 
to the necessity to deal with Lorentz indices. Let us first examine the expression in the second 
line of (5.5) with the four-point correlation function replaced by (5.7). According to (5.11), the 
Ad—functions corresponding to the i?—current and stress-energy tensor do not depend on Zi 2 
and admit the following representation MttJsJs ~ for 5 = 1 and S = 2, 

respectively. Then, we find from (5.7) 


/ OO 

dz^dz2 {T{l)T{2)Js{^)Jsm 

-OO 

/ OO 

dzf dz 2 F{z) 

-OO 

. (3|a;|4)2^ (4^4+)^ 


= {v^v^Y 


' —OO 
O \2r;^ 




(212) [dz^d^Y) [^12 ^(7)], 


(5.22) 


where in the second relation we applied (5.13). 

We observe the striking similarity of (5.22) with the first relation in (5.15). An important 
difference is however that the last relation involves the derivatives 113^4 and the matrix element 
(3|a;|4) defined in (3.24) and (5.12), respectively. They induce the dependence of (5.22) on the 
auxiliary spinors A* and Aj with i = 3,4. To obtain the correlation function in the second line 
of (5.5), we have to differentiate (5.22) with respect to A, and Aj. In this way, we arrive at the 
following differential operator 


Psidx3,dxY = (4 4 

xHa: 




iO'i'Ji 


(3)Ag‘(4)(9,3„. 


2=1 


iidx 




J3|a;|4) 


2S 


'3 ^4 


\S+2 


(5.23) 


where A“7(0 = It is easy to see that it does not depend 

on the auxiliary spinors and is given by a linear combination of powers of the differential operators 
dz 3 and dz^- Using (5.3) they can be converted into the differential operators dx^ and dx^- For 
5 = 1,2 the explicit expressions for Psidx 3 ,dxY are given below in (5.26). Notice that (5.23) 
only depends on points 3 and 4. 

Combining together (5.22) and (5.23) we find from (5.5) with the help of the identity (5.16) 

We observe that the only difference with (5.17) is the appearance of the differential operator 
Psidx 3 ,dxY that carries the dependence on the polarization tensor of the current. Upon the 
Fourier transform (5.1), this operator is replaced by Ps{q, —q) 

j 43:e“''Gs(a:;ni,n2) = (5-25) 


23 







where the function J^{x) is given by (5.19). Going through a lengthy calculation we find from 
(5.23) 


Ps=i{q, -q) = q^elei + ^(e“g“)(e^g^) = 2 q\^e, + 2{e>^q,){e,q '^), 

Ps^ 2 {q, -q) = 6(e:;^,:4^<?g)(e-::|Cgg) + 

= 24[(e^"<?^g.)(e^.g^<?") + l2q\^^Pep,q^q'^ - 


(5.26) 


where = ^e)l(cj^)“ and e^'' = 5 e“'^((T^)“((T^)|^ are the polarization tensors for currents of spin 
5 = 1 and 5 = 2, respectively. 

The very fact that the expression on the right-hand side of (5.25) factorizes into a product 
of the polynomial Ps{q,—q) and the universal scaling function P{x) immediately implies that 
the energy-energy correlations (5.1) normalized according to (5.20) cease to depend on the spin 
of the current and the polarization tensor. We therefore conclude that in AA = 4 superconformal 
theory the energy-energy correlations do not depend on the choice of the source 

EECj = EECt = EECo = EECl = . (5.27) 

In the A7 = 4 SYM theory, the function P{x) is known at next-to-leading order both at weak [21] 
and at strong coupling [18, 20, 22]. 


6 Conclusions 

In this paper, we have studied the four-point correlation functions of the conserved currents 
in AA = 4 superconformal theory. Contracting the Lorentz indices of the currents with auxiliary 
spinors we found that the correlation functions have a remarkably simple form - they are given by 
total spatial derivatives acting on some scalar functions, with the number of derivatives related to 
the spin of the currents. The scalar functions factor out into the product of a universal dynamical 
function and kinematical Ad—functions depending on the chiral auxiliary spinors and the space- 
time coordinates of the operators. We demonstrated that the requirement for the correlation 
functions to respect the conservation of the currents and to have correct conformal properties 
lead to powerful constraints on the Ad—functions. This allowed us to determine them for various 
correlation functions without doing any calculations. 

The obtained results for the Ad—functions revealed a surprising similarity with the known 
expressions for the scattering amplitudes in A/" = 4 SYM theory. Promoting the auxiliary spinors 
to new, twistor-like coordinates of the operators and assigning them a definite conformal weight, 
we found that the Ad—functions are built from very special blocks which have been previously 
identified as the simplest, four-point dual conformal invariants for the scattering amplitudes. We 
believe that this is not accidental and hints at the existence of an additional symmetry of the 
correlation functions. Our results are in agreement with the recent findings of [48] for a new class 
of A/" = 4 superconformal invariants. They admit a compact representation if expressed in terms 
of twistor variables analogous to (3.20) and their linear combinations describe the multiple-point 
correlation functions of A/" = 4 supercurrents in the chiral sector. It would be interesting to 
generalize such a representation to include the dependence on the anti-chiral variables 9. 
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It is natural to ask whether the higher-point correlation functions of the stress-energy tensor 
admit the same ‘derivative’ representation as the four-point functions. It is interesting to note that 
a similar representation exists for the two-point function (r(l)T(2)) {V^V2? [(12)V(x?2)'], 

but not for the three-point function since otherwise it would imply the vanishing of the two- 
point function, in virtue of the Ward identities [42], We recall that the distinguishing feature of 
the derivative representation is that the resulting expression for the correlation function satisfies 
the homogenous Ward identities (3.14). In general, the contact terms in the Ward identities 
are proportional to correlation functions with fewer points. Then, the absence of contact terms 
for the anomalous contribution to the four-point function is an immediate consequence of the 
protectedness of the two- and three-point functions in an A/" = 4 superconformal theory. Since 
the four-point function is not protected, we do not expect the higher-point correlation functions 
to obey homogenous Ward identities, thus making the derivative representation problematic. 

One of the byproducts of our analysis is the prediction for the four-point correlation function 
of stress-energy tensors in planar AA = 4 SYM theory at strong coupling. It is given by (3.2), (3.25) 
and (2.17) with the function <!>(«, v) replaced by its expression found in [37-39]. Via the AdS/CFT 
correspondence this correlation function is dual to the four-graviton scattering amplitude in AdSs. 
Due to the complexity of the corresponding Witten diagrams, such an amplitude has not been 
computed so far. The fact that the correlation function has the special form described above 
should simplify the problem and make the calculation feasible. 
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A Conventions 


A four-dimensional vector Xn can be represented by a 2 x 2 matrix 


Xaa — ^/i(<7^)oQ — 


X 


X X 


x'l = det\\Xaa\\= X~^x —XX. 


(A.l) 


We use the following conventions for rising/lowering indices 






aa _ ,a0 $a _ 


—X 

X 


where the completely antisymmetric tensors are normalized as 
satisfy the relations 


X 

—x~^ 

i2 


= ei 2 = e = = 


a/3^ _ cQ 

^ 7/9 ^7 ? 


^d0 . _ ra 

^ ^j0 - 


(A.2) 
1 and 

(A.3) 
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Defining the derivatives by {dx)°'°'Xaa = <5^, we find for Ao,d = (0, — 1) 


73/3 - "/ 3 ^/ 3 > 



Ao,d(5.)““ 



(A.4) 


It is straightforward to verify that under inversion = (x = x^ajx^ the derivatives 

transform as 



The generalization to an arbitrary number of derivatives is straightforward 


1 



(A.6) 


This relation can be proved by induction. Contracting the indices on both sides of (A.6) with 
e/3i/3i+i and using the identity we obtain 


/[□] = = (x2)3a;“i(5AiT'in)x5^ , /[tf] = . (A.7) 


where □ = 5^5^. 

B J\f = 4: superconformal generators 

In this appendix we elucidate the origin of the maximal number of derivatives recipe that we used 
in section 3.1 to compute various components of the super correlation function. 

It is well known that the special superconformal 5—transformations can be realized as a 
superposition of inversion and chiral super-Poincare Q—transformations 



(B.l) 


with the composite index A = {a, a'). The super-Poincare transformations generated by and 
Qa take the form 






(B.2) 


where x( = Xi and similarly for the other coordinates. The action of inversion on the 

supercoordinates looks as 




(B.3) 
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Combining these relations we obtain from (B.l) the global form of the S'—transformations 
x'l = x,(l - , y'l = (1 - + M ), 

, e'l = 9^ + x,(l - + i’yi) , (B.4) 


where x'l = Xi and we did not display the Lorentz and SU{A) indices for the sake of 

simplicity. One can verify that the generators of the transformations (B.4) are given by the 
operators S)? and S^ defined in (2.20). 

The expression for the supercorrelator (2.23) involves the differential operator which is 
related to the generator Sy^ by a similarity transformation (2.22). It is easy to see using (2.20) 
that 

^ (^2 2 2 2\2 
. [X-j^X2X^X^) 

~ (^2 2 2 2^^2^ Qb'0 I {^1^2^^^ 4 ) i (B-5) 

yxiX2X^x^) 

where in the second line we applied (B.l). It follows from the first relation that {S^^,^, Qa'a} = 0. 
Then, we apply (B.5) to get 




(^2^2^2„2\2 


T /^/4 7 "/ 2 2 2 2\2 

IQ /(X1X2X3X4) , 


(B.6) 


where S'^ = ^ Syj^ and similarly for Q'^. 

Let us now examine the action of on a test function f{x). Replacing the Q'—generator 
by its explicit expression (2.20) we find 


Q'V(x) = n [ E 

a' 


fix) 


1 

12 


E 

l<il,i2,«344<4 


ie^A^)^^^Ho^2^^Xno^^^i,)ZiOMa'c.Jix) 


a2(y.i . 


(B.7) 


where we used the notation for (0*9*)““ = and the SU{2) indices a and a' are 

lowered and raised in the same manner as the Lorentz indices, e.g, 0“ = e“^0b and 0a' = ^a'b'G^ 
(see Appendix A). Each term on the right-hand side of (B.7) involves four spatial derivatives and 
the same number of 0—variables. 

We can use (B.6) together with (B.7) to obtain an analogous expression for S"^f{x). To this 
end, we first apply inversion to both sides of (B.7) 

I Q'^f{x) = II {x\x\x\xA‘(t>ix) = {x\x^x\xA‘S'^4>ix) , (B.8) 

with (f){x) = I [(xfx|x|x|)^/(x)]. To evaluate II we apply the identity 

/[(09)“i“i ... (09)“-“"] = (x2)2(09)“i^i ... (09)“"^" x“) ... (B.9) 

[x ) 

that follows from (A.6) and (B.3). In this way, we arrive at 

= ^ E (45n)“^“'^(49*,)“^(0-*39*3)“3“2(0-*,9*J“J 

1 < 21 ,22,^3,24^4 

X (XiQaiaiixAa2»2(Xi3)aU^i4)al <f(x) • (B.IO) 
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Notice that the spatial derivatives in the first line do not commute with the product of x’s in 
the second line. Let us compare (B.IO) with the analogous relation in which we replace S' by its 
explicit expression (2.20) 


?'V(x) = n[ E 


a'.a 1 < 2<4 


+ 0<9^M^of^ 


(j){x) 


(B.ll) 


We observe that in order to reproduce (B.IO) it is sufficient to neglect the second term inside the 
brackets in (B.ll), and then move all spatial derivatives to the left of all x—dependent factors. 
It is this shortcut that we used in section 3.1. 
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